February 6, 2008 22:41 Advances in Physics perspective 



Advances in Physics 

Vol. 00, No. 00, Month-Month 200x, 1-25 



o 
o 

^ . 

O ■ What makes a crystal supersolid ? 

Nikolay Prokof 'ev, Department of Physics, University of Massachusetts, Amherst, 
C^ . MA 01003, USA 

^ I (Received 00 Month 200x; In final form 00 Month 200x) 

Wpor nearly half a century the supersolid phase of matter has remained mysterious, not only eluding 
Gsxperimental observation, but also generating a great deal of controversy among theorists. Recent 
Ldiscovery of what is interpreted as a non-classical moment of inertia at low temperature in solid 

. ^^l\e [1-4] has elicited much excitement as a possible first observation of a supersolid phase. In the 

' (^wo years following the discovery, however, more puzzles than answers have been provided to the 
^ fandamental issue of whether the supersolid phase exists, in helium or any other naturally occurring 
jrtpondensed matter system. Presently, there is no established theoretical framework to understand the 
*-Hbody of experimental data on ''He. Different microscopic mechanisms that have been suggested to 
^underlie superfluidity in a perfect quantum crystal do not seem viable for ''He, for which a wealth 
I o f experimental and theoretical evidence points to an insulating crystalline ground state. This per- 

' Ofepective addresses some of the outstanding problems with the interpretation of recent experimental 
^Observations of the apparent superfluid response in 'He (seen now by several groups [5-7]) and 
Qdiscusses various scenarios alternative to the homogeneous supersolid phase, such as superfluidity 
Q_Jinduced by extended defects of the crystalline structure which include grain boundaries, dislocations, 

I— Anisotropic stresses, etc. Can a metastable superfluid "glassy" phase exist, and can it be relevant to 
some of the experimental observations ? One of the most interesting and unsolved fundamental ques- 

^— Htions is what interatomic potentials, given the freedom to design one, can support an ideal supersolid 
^^hase in continuous space, and can thoy be found in Nature. 

OS ■ 

o\ : 

Introduction 

^iThe textbook notion of a perfect crystal at T = 0, is that of a periodic array of 
"^unit cells, all comprising the same integer number of particles, = Jq drp{r), 

Svh ere is the unit cell volume and p{r) is the average particle density pro- 

Qile. For simplicity, let us consider a single component crystal similar to ^He. 
'^We assume that the elementary constituents (atoms or molecules) can be re- 

Cngarded as structureless particles; furthermore, we assume that they obey Bose 

^statistics. 

. The supersolid phase (SFS) can be generally defined as one that combines 
. S^rystalline properties, such as shear modulus and broken translation symmetry, 
^^ith frictionless mass transport through the solid bulk. The striking, simulta- 
j^eous presence of solid and superfluid properties in the same condensed matter 
■ ' system, will result in a number of phenomena that defy our everyday experi- 
ence. One such phenomenon is schematically depicted in Fig. 1, showing a SFS 
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Figure 1. Supersolid crystal (represented by periodic arrays of dots) in coexistence with the 
superfiuid liquid of lower density (shaded regions on top of the U-shaped closed tube) in the 
gravity field. If initial conditions are such that supersolid levels to the left and to the right are at 
different height, the system will start oscillating by melting on one side and recrystallizing released 

density on another side with coherent dissipationless mass flow through the otherwise static 
crystalline structure. [The assumption is that the flow velocity remains below the critical value, or, 
h is small enough, and the supersolid-liquid interface is rough.] 

sample (region with periodic arrays of dots) placed inside a U-shaped vessel, 
in coexistence with the superfiuid liquid. At T = undamped oscillations of 
the solid levels in the two sections of the U-shape vessel will be observed, as 
long as the maximum flow velocity remains below the critical value and the 
liquid-solid interface is rough. 

Another landmark of SFS behavior (see, for instance, Ref. [8]), is the reduc- 
tion of the moment of inertia of a solid sample with respect to its classical 
value (commonly referred to as Non-Classical Rotational Inertia, or NCRI) 

/(T) =/ciass(l-ris(T)/n) . (1) 

Here ns{T) is the superfiuid density, and n is the particle density in the system. 
The NCRI effect can be observed by enclosing a known amount of solid in a 
vessel, which is then set in rotation about its axis. The moment of inertia 
of the system is related to the (measurable) resonant period of oscillation of 
the rotating system. A drop of the moment of inertia at low temperature, is 
interpreted as decoupling of part of the solid (the superfiuid fraction) from the 
rotation. 

The history of ideas on how superfluidity can occur in a crystalline solid is 
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Figure 2. A sketch of the ^Hc phase diagram. At low temperature the heUum superfluid undergoes 
a weak first-order phase transition to the hexagonal close-packed (hep) solid phase at pressure 

Pm « 25.6 bar. 



quite old. On very general grounds, one expects that the "textbook" crystal 
ought to be insulating, i.e., with vanishing superfluid density Ug, and off- 
diagonal correlations decaying exponentially in space. Penrose and Onsager 
(PO) were the first to argue in favor of this point of view [9, 10], starting 
from the picture of atoms localized around their equilibrium lattice positions. 
However, their argument has never been regarded as the final word in the 
discussion of whether the superfluid solid phase of matter is possible. C. N. 
Yang made a comment that long-range phase correlations in the solid state 
could occur in systems characterized by a high degree of atomic delocalization 
[11]. Helium is unique in this regard; quantum properties of light helium atoms 
in combination with the relatively weak interatomic potential (the well depth 
is only ~ 10 K) prevent them from making a crystal at zero pressure, see 
Fig. 2. An objection to the PO treatment was that regarding atoms as being 
localized around lattice points underestimates the role of exchange processes, 
which may allow for the (perhaps remote) possibility of superflow, even with 
an integer u [8, 12]. 

Another proposal for the supersolid phase has its origin in the observa- 
tion that there is no fundamental reason why crystals should necessarily be 
commensurate, i.e., feature an integer number of particles per unit cell, on 
average. A dilute gas of vacancies, or interstitials, which occur in all solids at 
flnite temperature and ought to be highly mobile in a helium crystal, could 
be also present in the ground state giving rise to Bose-Einstein condensation 
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and superfluidity at low temperature [13, 14]. 

Though certainly plausible, most theoretical scenarios leading to the SFS 
phase are phenomenological and qualitative when it comes to making predic- 
tions for a particular Hamiltonian; they have never been supported by robust 
microscopic calculations, including all quantum-mechanical effects and based 
on realistic models of crystals. ^ Still, their appeal is so strong that they in- 
spired decades of active research and continue to do so. 

Remarkable manifestations of quantum behavior on a macroscopic scale in 
supersolids are expected to render the observation of the SFS state rather 
unambiguous, as long as it occurs under conditions (e.g., temperature) acces- 
sible to the experiment. Such an observation has eluded experiments for some 
thirty years (for a review of the experimental evidence prior to 2004, see, for 
instance, Ref. [17]); it was only in 2004 that the first convincing evidence of 
a SFS phase of ^He was published. The phenomenon was first reported for a 
solid embedded in a porous Vycor matrix by [1] Kim and Chan (KC); shortly 
thereafter, it was extended to the system embedded in a matrix of porous gold 
(with a characteristic pore size two orders of magnitude greater than in Vycor), 
and eventually to the bulk crystal [2]. In all cases, the observed NCRI at low 
temperature was identified as the onset of supersolidity in ^He. It should be 
emphasized that the supersolid interpretation of the data is supported by the 
velocity /amplitude dependence of the NCRI signal and the crucial test that 
the NCRI effect goes away (is suppressed by two orders of magnitude) when 
the flow channel around the rotation axis is blocked. These results constitute 
a major problem for the non-superfluid kinetic interpretation. 

The Kim and Chan discovery sparked a renewed effort in the investigation of 
the SFS phase of matter. At the time of this writing, it seems fair to state that 
things have turned out to be considerably more complicated (and thus more 
interesting) than expected. As the initial experiment by Kim and Chan is being 
repeated by other groups, and concurrently different experiments and reliable 
microscopic (first-principles) calculations, based on state-of-the-art computa- 
tional techniques are being carried out, consensus is building that the early 
microscopic scenarios of supersolidity, described above, are not viable to ex- 
plain the body of experimental data. For example, there is now strong theoret- 
ical and experimental evidence that ^He is a commensurate crystal, and such 
crystals made of a single species of particles obeying Bose statistics (e.g., ^He 
atoms) are always insulating (more precisely, the commensurate SFS phase 



^ A major current area of investigation, in tlie broader context of the study of exotic phases of matter, 
focuses on quantum lattice models. Supersolid phases liave been predicted to occur in some of tliese 
models [15, 16] and may soon be observed experimentally, e.g. in optical lattices. Albeit of clear 
fundamental interest, this topic is not part of the discussion in this paper. We also do not discuss 
multicomponcnt cases in which one component forms a solid and, thus, acts as an external periodic 
potential for the other. Such cases are expected to behave similarly to lattice models. 
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has zero probability of being observed) . Measurements and calculations of va- 
cancy and interstitial properties indicate that the Andreev-Lifshits-Chester 
scenario of vacancy-induced supersolidity does not occur in solid helium, even 
out of equilibrium. Therefore, attention is now shifting to a variety of inhomo- 
geneous scenarios, that may account for many puzzling (and often apparently 
conflicting) experimental results. 

This paper is organized as follows: we shall first review the theoretical frame- 
work for supersolidity, including the early suggestions as well as some more 
recent proposals of microscopic mechanisms. They will be discussed in the light 
of both early and recent analytical and numerical results. We shall then ex- 
amine the core of experimental data, which at this time includes observations 
of NCRI as well as flow experiments [18] and the most recent experiment on 
grain boundary superfluidity [19], and propose that, while a definitive theoret- 
ical explanation of the observed phenomena is not yet available, there are good 
reasons to look at disordered and/or inhomogeneous scenarios. Specifically, we 
shall illustrate how the effect of disorder can hold the key to interpreting an 
important part of the phenomenology. Finally, we shall outline some possible 
directions of experimental and theoretical investigation. 

This perspective is a quick response to the fast developing topic, not a com- 
prehensive review. I apologize if some important contributions to the field were 
overlooked. Given a broad range of opinions and opposing points of view on 
the subject, there is certainly a degree of personal bias towards which argu- 
ments and results are more convincing and should be given more "weight" 
in the discussion. The hope is that the story is provocative enough to stimu- 
late additional research efforts to understand the remarkable phenomenon of 
supersolidity. 



2 Supersolid ground state 

Often the discussion starts with the proposal for the supersolid ground state 
wave function though it is certainly allowed by the theory that in a par- 
ticular material SFS exists only as a finite-temperature phase. In the Penrose- 
Onsager picture of atoms localized around equilibrium points {Ri}, the ground 
state wave function is simply a product of localized Wannier orbitals 




p j=i 
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where the sum is over all permutations P of particle labels. Using the 
definition of the single-particle density matrix at zero-temperature 

n(r,r')= J . . . J dr2drN'^G{r,r2, ■ ■ ■ ,rN) ^cir' ,r2, . . . ,rN) , (3) 

and simple properties of the variational state one readily obtains (by 

substituting ^^f^^ for '^a above) 

n{r,r') = ^Y^cp{R,-r)cp{R,-r') . (4) 

i 

Since in the |r — r'j oo limit either r or r' has to be far from the equilibrium 
lattice point Ri, the density matrix decays to zero exponentially fast for lo- 
calized Wannier orbitals. This immediately implies no off-diagonal long-range 
order (ODLRO) and thus no superfluidity in the system. 

Here, we adopt the following definition of ODLRO: we say that it is present 
in the system, if the integral 

J n(r, r')dr oo , (Def. ODLRO) (5) 

diverges in the thermodynamic limit. This definition is different and more 
general than that based on finite condensate density (no), which requires that 
the integral be proportional to the system volume, in the thermodynamic 
limit [8,10,20,21]. It makes perfect sense to define ODLRO in such a way that 
superfluidity automatically leads to ODLRO and vice versa. In particular, 
Eq. (5) allows one to talk about ODLRO in two dimensional (2D) systems at 
finite temperature (often referred to as topological order [22]) for which the 
condensate density is zero. 

In this article, we restrict ourselves to the single-particle superfiuidity. In 
general, one has to look at the many-body density matrix to see if there is 
ODLRO in some finite-order channel. For example, it is certainly possible 
that M atoms form M-molecules and superfluidity happens than at the molec- 
ular level — in this case one has to study the M-particle density matrix. Or, 
molecular hydrogen would be a good example here, by neglecting the internal 
structure of molecules and treating them as point particles (probably at the 
expense of introducing complex effective many-body interactions) ODLRO can 
again be discussed in the single-particle "molecular" density matrix. There is 
no fundamental difference between the two pictures. 

The superfluid density (formally a tensor) is a linear response coefficient 
which controls the free-energy increase in a system with periodic boundary 
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conditions (BC) when these BC are twisted, i.e. if Lq, is the hnear system size 
in direction a = x,y, z we require that ^'(. . . , Va+L^, ■ ■ ■) = e"^" • • , ra, ■ ■ •) 
(there are other definitions of Ug, e.g. through the NCRI, Eq. (1), which can 
be shown to be identical to the one given below) 

aa' 

Here V = L„ is the system volume. Twisted BC can be avoided by intro- 
ducing instead a gauge phase with the gradient Vi;^ and writing the free-energy 
density (on the largest scales) as 



/(0)-/(o) = E^^(V'^)«w)" 

aa' 



2m 



Using "connectivity" arguments [8, 20, 21] one can show that the system is 
sensitive to twisted BC only if it has ODLRO defined in Eq. (5). 

The wave function (2) takes into account zero-point motion and local ex- 
change processes between identical particles (these effects are included in the 
shape of the Wannier orbitals); nevertheless, it describes an insulating state, as 
long as ^p{Hi — r) are localized. By moving away from two major assumptions 
made in the construction of the insulating ground state (2), one arrives at 
two, apparently quite different possibilities for the supersolid wave function. If 
the role of exchange is underestimated in localizing particles around different 
lattice points, then a better variational state might be of the BEC form 



with macroscopic condensate in the single-particle state ^jLi fi^j ~ 

The BEC state (as well as '^q''^^) can be further improved by multi- 
plying it by a many-body correlation factor of the Jastrow form, namely 
exp{— u{ri — rj)}, whose purpose is to account for the strong repulsive 

core of the interatomic potential, and suppress the condensate fraction (no) 
to a value much smaller than unity (no is already below 10% in the superfiuid 
^He at zero pressure [23-25]). By construction, the Jastrow factor ought not 
alter, nor eliminate the underlying lattice structure. 

Equation (7) offers an alternative view of the supersolid state, namely that 
of a superfiuid phase with a density wave modulation. However, it suffers from 
a fundamental shortcoming, in that the number of lattice points Nl is arbi- 
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trarily assumed to be the same as the number of particles. In a superfluid the 
density-wave parameters of any variational ansatz are independent thermody- 
namic variables, to be selected through energy minimization. At this point, 
we recall that superfluidity and ODLRO imply that part of the system matter 
is characterized by a classical field component, for which the notion of parti- 
cle number can be safely ignored without loss of generality. [This component 
represents low-momentum states with large occTipation numbers.] Thus, the 
condition = N x integer can be satisfied only by accident, which is to 
say that a commensurate supersolid in continuous space has a probability of 
zero measure to be found in Nature. As a result, Eq. (7) can not describe 
the generic ground state of a realistic system [26]. In order to illustrate more 
vividly how ODLRO relates to the presence of gapless vacancies and/or in- 
terstitials and incommensurability of a single component solid in continuous 
space, consider the simple example of a large-amplitude (classical) standing 
wave of the electro-magnetic field in the typical optical table experiment. In 
this case, the question would never be asked whether the number of photons 
per wave period is an integer. 

Another possible modification of Eq. (2), one that retains the picture of 
localized orbitals, consists of allowing the number of particles N to be less 
than the number of lattice sites A''^;,, i.e., assuming that the ground state has 
no energy gap for the creation of vacancies^ 



where the the sum is over all possible sets of N lattice points out of Nl 
available. Now, the outcome of the calculation for the density matrix is quite 
different because in the sum over i,j 



one can always find terms with Ri and Rj close, respectively, to r and r' , no 
matter how large their separation. If, in the spirit of a tight-binding approach, 
we move from a description in terms of continuous space variable r and r', 
to that of lattice site positions for the particle, by introducing site creation 
(and, analogously, annihilation) operators bj = J dr (p{Ri — r) bt, then Eq. (9) 



-•^Thc same argument can be used for intcrstitials, although typically their energy cost is higher than 
for vacancies; henceforth, we shall implicitly make this assumption consistently throughout the paper 
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can be recast in the following form (taking lattice translation invariance into 
account) 

n{l) = n{k, k + l) = Y,{hlhk+i) = ^^^^ ^ . (10) 
k k 

The final result is consistent with the picture of BEC of a non-interacting 
vacancy gas, with dimensionless lattice concentration n^, introduced by An- 
dreev and Lifshitz [13]. One has to assume next that the effective interaction 
between gapless vacancies is repulsive, in order to ensure system stability (i.e., 
to prevent the vacancy gas from collapsing). 

As was noted by Reatto [27] and Chester [14], it is not at all necessary to 
break translational invariance explicitly in the SFS state with vacancies, by 
specifying the set of equilibrium particle positions {-Rj}. An analogy between 
the Jastrow wave function 



{J) 

Q oc exp 



N 
i<j=l 



u{ri - Vj) 



(11) 



and the partition function of a classical system of interacting particles with 
pairwise potential v^'^\r) = Tu{r) at finite temperature, suggests that among 
(11) there arc variational states which feature spontaneous crystalline order, 
as classical systems are known to crystallize at sufficiently low temperature. 
States (11) are also superfluid since they feature ODLRO, and incommensurate 
[27] . Regarding the last property, one may observe that, due to lack of vacancy- 
interstitial symmetry in continuous solids, the probability of having = Nl 
is of zero measure [26] (an obvious statement for a classical solid at finite T). 

Based on all of the above considerations, one might be led to thinking that 
the Jastrow wave function may hold the key to the microscopic understanding 
of supersolidity. However, there are fundamental problems associated with 
Eq. (11), chiefly that it can not possibly describe realistic solids of systems with 
short-range interactions (such as helium), for which the condensate fraction is 
tiny (less than 10~^). The reasoning supporting this criticism is offered in the 
next Section. 

The most important conclusion of this Section is that only zero-point va- 
cancies or intcrstitials can make an ideal crystal supersolid. There is no other 
ideal-crystal scenario for the SFS state though microscopic mechanisms lead- 
ing to the incommensurate crystalline groundstate might be very involved. 
However, the story does not end here, since we did not question yet what 
happens when quantum crystals are not perfect, i.e. contain defects such as 
dislocations and grain boundaries, or even loose their crystalline order com- 
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pletely and form an amorphous solid. It appears, that topological defects and 
glasses offer an alternative approach to super solid phenomenon. In what fol- 
lows we examine which of the two possibilities is most likely to operate in 
helium. 



3 Insulating hep crystals of "^He. 

Existing theories can not answer the question at what pressures and tempera- 
tures one has to look for the SFS phase, or even whether the SFS phase exists 
at all in a given Hamiltonian, e.g. in helium. Several aspects have emerged as 
necessary conditions that the SFS state has to satisfy, in order to exist; conse- 
quently, any of them can be used to probe hep crystals of helium theoretically. 
[At the time of this writing, there is no experimental evidence of structural 
transitions in solid helium at low temperatures T < IK and moderate pres- 
sures, say P < 200 bar.] By definition, the SFS state has a finite superfluid 
density n^. Since superfluidity and ODLRO (as defined above) must be present 
simultaneously, one can also look at asymptotic properties of the single-particle 
density matrix, which is nothing by the zero time limit n(r) = G(r, —0) of the 
Matsubara Green function 

G(r , r)=V~'l dr'{ T{h,,^Ar) 4' (0)} ) • (12) 

Here (...) stands for the thermal expectation value, T is the time-ordering 

operator, —(3/2 < r < 13/2, and h\.{T) is the Bosc particle creation operator in 
the Matsubara representation. For r < (r > 0) one is computing the Green 
function for a vacancy (interstitial atom). If, over large distances and/or long 
times, n(r) and G(r, r) decay exponentially, then the state is non-superfluid, 
i.e., insulating. Moreover, if n(r) and G(r,r) are temperature-independent up 
to a certain length/time scale, then one can claim that ground state properties 
on the corresponding scales are being explored (this is most readily seen in 
the path-integral framework) . Finally, by performing simulations in the grand 
canonical ensemble one can investigate whether the solid is commensurate or 
has gapless vacancies/interstitials in the ground state. 

All of the above criteria have been used in recent first-principles simulations 
of hep helium crystals. They are based on the standard pairwise interatomic 
potential for helium [28] which essentially did not change in the last decades 
and is known to capture all important properties of condensed and gaseous 
helium with the relative accuracy of the order of one percent. The study of 
exchange cycles [29] reveals that they are extremely rare in the solid phase 
and indicative of the insulating behavior [30] . One may recall that macroscopic 
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Figure 3. Simulation results for the ideal hep crystal of A'^=800 atoms at low temperature 
T = 0.2 K. Left panel [32]: the density matix close to the melting curve, P = 32 bar (n = 0.0292 
A"'' ), and at elevated pressure P = 155 bar (n = 0.0359 A"'' ). The solid line is representing an 

exponential decay. Right panel [38]: the zero-momentum Green function at the melting density 
no =0.0287 A~^ . Symbols refer to numerical data, solid lines are fits to the long-time exponential 
decay. The given numerical values are the interstitial (Aj = 22.8 it 0.7 K) and the vacancy 
(Av = 13.0 it 0.5 K) activation energies, inferred from the slopes of G. The inset shows the 
vacancy-interstitial gap -Egap = Ai + Av for different system sizes. 



exchange cycles are necessary for superfluidity, since 

mT{W'^) 

n, = 

dL 

where m is the helium mass, d is the dimensionality, L is the linear system 
size (periodic boundary conditions are assumed), and W = {Wx, Wy, Wz) are 
winding numbers, counting how many times exchange cycles involving many 
particles wind around the system boundaries [31] . To simplify, we ignore here 
the tensor structure of Us in the crystal. Initial concerns about ergodicity 
problems of simulations probing large exchange cycles [30] were eliminated 
in subsequent path-integral Monte Carlo (PIIVIC) simulations based on the 
Worm Algorithm [32], which established that the superfluid fraction n^/n is 
unobservably small (below 10~^), even in finite-size crystals comprising only 
800 atoms. 

More detailed evidence for the insulating behavior of helium crystals is 
provided by computations of the single-particle density matrix and the zero- 
momentum Matsubara Green function [32,38]. In Fig. 3, we clearly see an 
exponential decay of the one-body density matrix n(r) at large distances, 
both in the vicinity of the melting curve (pressure P = 32 bar) and at the 
pressure of P = 155 bar. The observed behavior is independent of temperature 
for T < 1.5 K, and there is good agreement between independent calculations 
of n{r) performed by different groups at the melting pressure [33]. 

From the exponential decay of the zero-momentum G(k = 0, r) (left panel 
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of Fig. 3), one can obtain an activation energy A = Ay ~ 13 K for a vacancy 
(Ai « 23 K for an interstitial) by studying 



for r < (r > 0) . We note that these energies are quite large compared to the 
characteristic temperature of 0.2 K of the KC experiment, all but ruling out the 
possibility that thermally excited vacancies may be behind the supersolid phe- 
nomenon (it should also be noted that Ay and Ai increase with pressure) [38] . 
This conclusion is consistent with the absence of thermally activated vacan- 
cies at low T, as inferred from impurity mobility experiments (summarized 
in Ref. [34]), as well as high-precision studies of the liquid-solid phase dia- 
gram [35] . One may note that simulations of gaps in a system comprising 800 
atoms are essentially probing the thermodynamic limit, as finite-size effects 
are negligible (see inset of left panel of Fig. 3). 

The analysis of the X-ray data at low temperature in terms of activated 
vacancies [34] was questioned in Ref. [40]. Specifically, the formula 



with reported activation energies, in combination with the calculated [36] va- 
cancy effective mass, can not account for the observed change of the lattice 
constant (the effect would be too small). However, since energies of optical 
and acoustic phonons at the Brillouin zone boundary are also ~ 10 K, one 
can not exclude a conventional explanation in terms of crystal anharmonicity. 
It seems, that all one can state with relative confidence is that the vacancy 
concentration is below 1%. 

In the most conservative approach, the above-mentioned theoretical results 
show that commensurate hep crystals of ^He are insulating and locally stable 
since single vacancies have a finite energy cost. This does not exclude, strictly 
speaking, a remote possibility that the true ground state may have a finite, 
albeit small concentration of vacancies [39,40]. Such scenarios, however, are 
ruled out by strong attractive interaction between vacancies; indeed, when 
vacancies are forcibly introduced in the crystal at any finite concentration (in 
a "computer experiment" ) , they immediately phase separate from the crystal 
bulk [38]. Under realistic experimental conditions, they will quickly anneal 
at grain boundaries and dislocations; if the latter are absent, then vacancies 
will form liquid droplets (at the melting curve) or dislocation loops. In any 
case, the net result is that the commensurate insulating hep crystal is the true 
ground state. 

In addition to the original picture of the dilute zero-point vacancy gas 



G(k = 0,r) cx e 



r| — ^ GO, 
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[13, 14], over the past few years there have been several theoretical proposals 
for the ideal supersolid phase of helium [39,41-44], differing in the microscopic 
mechanisms of supersolidity. Some of the early works suggesting commensu- 
rate SFS phases, were subsequently amended, to respect the theorem that such 
phases should necessarily feature gapless vacancies (not necessarily in the di- 
lute gas phase). In this respect, Refs. [41,43] are still in error. An interesting 
idea was proposed in Ref. [45], that hep crystals may become SFS only un- 
der anisotropic stresses, which are likely to occur in the experiment. However, 
simulation results for hep crystals at the melting curve, show that vacancy and 
interstitial gaps hardly change, even under very large anisotropic stress [46]. 
It appears at this point that phcnomenological scenarios and mean-field theo- 
ries of supersolidity do not work for "^Hc, and one has to seek the explanation 
of recent experimental results outside the paradigm of the homogeneous hep 
crystal. 

Before reviewing crucial experimental facts, let us discuss the most impor- 
tant shortcomings (including some unphysical properties) of Jastrow, as well 
as of the related "shadow" [47] variational wave functions, in the context of 
helium. Calculations based on these wave functions predict both (i) a very 
small, but finite, condensate fraction for solid ^Hc near the melting curve 
(below 10"*^) and (ii) large activation energies for vacancies in solid helium 
samples comprising several hundred atoms [48] . In the presence of short-range 
interactions among atoms, the physics of vacancy formation is local in nature. 
All of the energy contributions quickly approach their thermodynamic limit 
value, as the system size is increased. Indeed, the kinetic energy of localiza- 
tion in a volume is proportional to (this contribution is not present 
in a system with periodic boundary conditions). The direct coupling of the 
removed atom with the rest of the system converges as for the Van der 
Waals interaction; the same law describes convergence of the deformation en- 
ergy. The fast disappearance of finite-size effects as L is increased, is observed 
in essentially all numerical simulations of condensed helium (see, e.g. inset in 
Fig. 3), including those based on the variational approach. On the other hand, 
vacancies are gapless in a state described by a Jastrow wave function by con- 
struction! It means that if variational calculations were taken to larger system 
sizes they would have predicted negative corrections to the vacancy energy 
which inerease, not decrease, with L, making A/ negative for 77,0-^'^ > 1- This 
unphysical behavior is only possible if irreducible long-range multi-particle in- 
teractions, which can not be reduced to the effective chemical potential shift, 
are present in the system — something that is not observed in Nature. 

The other problem with the Jastrow variational ansatz (and, of other such 
wave functions as well) is that it does not easily describe the physics of phase 
separation. Because, as shown above, a dilute Bose gas of vacancies in helium 
has attractive interactions and thus is an unstable, phase separating system. 
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any discussion of BEC and superfluidity of the uniform vacancy gas under 
such circumstances is meaningless. 

4 ^He "ice cream" 

In this Section, we shall argue that many experimental facts suggest that the 
observed supcrsolid phenomenon in ^Hc is induced by defects, or disorder, 
in the torsion oscillator samples. This point of view is supported by several 
numerical simulations, as well as by the recent direct observation of the grain 
boundary superflow [19]. 

Let us assume that the observed drop in the resonant period of the tor- 
sion oscillator can be attributed to a non-zero superfluid response at a given 
frequency. This seems plausible, as the period drop or, equivalently, mass de- 
coupling increases as the oscillator amplitude Aosc and velocity I'osc decrease. 
[To understand why decreases with Vosc, one has to consider the non-linear 
response of vortex loops and pinned vortex lines to the flow.] In normal vis- 
cous media, an opposite correlation between the mass decoupling and Vosc is 
expected. One may also attempt to understand the NCRI effect in terms of 
kinetic relaxation or mass redistribution in the cell [49,50]. However, in helium 
samples, these explanations face the crucial blocked-channel test [49] -only su- 
perflow can be sensitive to the sample topology on macroscopic length-scales. 
It is also hard to believe that mass redistribution kinetics in solid helium is 
reversible at low temperature, i.e. the signal would not be reproducible upon 
heating and cooling (assuming that T < 1 K in the cycle to avoid possible 
annealing effects). 

With the superfluid response in mind, we observe that many features in the 
data simply do not fit into the homogeneous supersolid picture. In particular: 

Dependence of the superfluid density on temperature 

It is a well established fact that in a continuous normal-sup erfluid (N-SF) tran- 
sition, the superfluid density dependence on temperature near the transition 
point at Tc is given by (T — TcY, where u « 0.671(1) is the correlation length 
exponent for the 3D XY-universality class [51]. So far, no exception to this 
law, which predicts an infinite derivative dug/dT at has been found, either 
theoretically or experimentally, and the same is expected for the transition to 
the SFS state [45]. On the other hand, in all torsional oscillator experiments 
displaying mass decoupling at low T, the rise of n<j(T) at low T starts off 
with zero derivative near Tc, in turn rendering a precise determination of Tc 
rather ambiguous. The most straightforward explanation for such a gradual 
increase of the superfluid density would be a broad distribution of local transi- 
tion temperatures T'c(r) within a strongly inhomogeneous sample. The shape 
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of the Tig (T) curve is then controlled by the probability distribution of transi- 
tion temperatures. Diffraction experiments might provide information on the 
crystal quality in experimental samples but these were not reported yet. 

Specific heat anomaly 

The other puzzling feature, which does not conform to the established pic- 
ture of the continuous N-SF transition, is the absence of the specific heat 
maximum. Given experimental uncertainties [3], the possible amplitude of the 
specific heat anomaly is orders of magnitude too small for the observed amount 
of the superfluid density [52]. Moreover, it was found that the dependence of 
the specific heat on T is linear in the supersolid regime. A linear Cy{T) 
may originate from one-dimensional Luttinger-liquid type structures, e.g., su- 
perfluid dislocation cores [52] and ridges between grain boundaries [26,46], or 
just from a collection of two- level systems in the amorphous sample [53]. 

Dependence on geometry 

If the supersolid phenomenon observed in recent experiments is indeed a ho- 
mogeneous bulk effect, it is difficult to explain the observed dependence of the 
microscopic^ parameter Ug on the geometry of the experimental setup . The 
scatter of the reported values of Ug (spanning almost a decade) , depending on 
annulus or open geometry in the KG experiment, can be in principle explained 
if the sample quality, and thus the distribution of Tc, is strongly dependent on 
the cell geometry (certainly a plausible argument). 

Dependence on sample "history" 

Rittner and Reppy [5] have reported the "elimination of supersolid by an- 
nealing" , in an experimental setup very close to that of Kim and Ghan. The 
other important observation, is that the torsional oscillator quality factor, or 
inverse dissipation rate, is higher in normal samples. This result should be re- 
garded as direct evidence in favor of disorder-induced supersolidity. However, 
such annealing effects have not been confirmed by Kim and Ghan [4] , and the 
explanation for this discrepancy is lacking, at the time of this writing. 

^He effect 

Bosonic superfluidity is a very robust phenomenon, and one does not expect 
any significant changes in the superfiuid properties when ^He impurities are 
added, at a concentration ns << n^. It is well established [54-56] that ^He sub- 
stitution atom in the ^He solid matrix is described by a tight-binding model, 
with a tiny hopping amplitude J ^ 10~^ K. The small value of J is consistent 



^Microscopic is defined here as occurring on length scales much larger than the correlation length 
but much smaller than the system size 
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with the characteristic energy for nuclear magnetism in sohd ^He, which is 
beUeved to be due to exchange of heUum atoms [57]. This result is important 
in several ways. First, it proves that ^He atoms do not induce vacancy for- 
mation in their vicinity since the ^He-vacancy complex would be extremely 
mobile (the vacancy hopping amplitude is nearly four orders of magnitude 
higher [36]). Second, a small value of J makes ^He atoms extremely sensitive 
to virtually any type of crystalline disorder, since defects will certainly create 
tight bound states for ^He. Even a weak deformation potential between ^He 
impurities is sufficient to bind them [54,58], or localize dilute ^He solutions [59]. 

It is therefore very surprising, that both the superfiuid fraction and the 
onset of the superfiuid response observed in KC experiments, are significantly 
altered when ^He impurities are added, even if the concentration is as small 
as few ppm, (or even few ppb !). The mismatch between Ug and ns is so large 
that it casts serious doubts on the homogeneous supersolid scenario [26] , even 
on phenomenological grounds [60]. 

In a disordered sample, the lighter ^He atom will try to minimize its kinetic 
energy by binding to static defects which have lower local particle number 
density, such as dislocation cores, grain boundaries and ridges between them 
(vacancies are mobile and it is not clear whether they can form a bound state 
with ^He atoms). One thus expects an accumulation of ^He on defects at low 
temperature; it is plausible that the local concentration of isotopic impurities 
may be fairly high in certain regions of space. The other, somewhat speculative, 
effect of '^He substitution might be that the quality of samples, i.e. the amount 
of disorder, is itself a function of ns. For example, by concentrating at the 
perimeter of the microcrystal, ^He atoms may inhibit fast crystal growth, 
thereby helping to make the "ice cream" . 

If we now couple these considerations with the theoretical conclusion that 
only disorder can be responsible for the supersolid phenomenon in helium, we 
find a very subtle interplay between the superfiuid response and ^He content. 
Imagine for a moment, that superfluidity occurs along grain boundaries, while 
■^He is concentrated at the boundary ridges. Assuming that '^He atoms suppress 
superfluidity, we arrive at the picture where ^He-rich ridges form a Josephson 
junction network completely surrounding grain boundaries. In this scenario it 
is conceivable that even a tiny amount of ^He (per volume) might be relevant 
for supersolidity [26]. 

Dependence on pressure 

The obvious expectation is that supersolid properties ought to be suppressed, 
as pressure is applied to the sample. Applying pressure indeed drives down 
the superfiuid transition temperature Tc. Practically all known properties of 
helium crystals are consistent with the overall tendency of helium atoms to be- 
have more classically at higher densities (see, e.g. Fig. 3). Surprisingly, the KC 
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data [4] show hardly any pressure dependence for Ug and Tc, up to P ~ 65 bar-^; 
at higher pressures, the supersohd signal weakens, and apparently disappears 
at P > 160 bar. At this point one may wonder whether experimental data on 
pressure dependence should be attributed to properties of the same sample, 
as if the pressure was changed by squeezing the volume, or different point 
represents properties of different samples. If supersolidity is directly related to 
the sample quality, which, in turn, depends on initial conditions and the solid- 
ification protocol, then data become very ambiguous, and their interpretation 
intricate. 

While the above mentioned experimental facts do not fit into the homoge- 
neous supcrsolid crystal picture, there are others which are hard to reconcile 
with any existing theoretical framework. For example: 

"Critical velocity" 

Superfluid decoupling in the torsion bob has the most unusual dependence 
on the oscillator amplitude/velocity. Though in the current setup velocity 
and amplitude are strictly related, it is assumed that forces equilibrate in 
the sample fast enough to attribute all changes in the supcrflow to the ve- 
locity V. It is found that Ug saturates at macroscopically small velocities, 
^^sat ~ 10 fim/s ~ h/mL (related to just few circulation quanta). These 
values clearly have nothing to do with the naive sound velocity estimate 
c ~ ^UgU jm ~ 10 m/s based on typical helium parameters. The discrep- 
ancy is so large, that there is virtually no room left for explaining v^^x. using 
microscopic mechanisms. On the other hand, if extended defects in the super- 
fluid phase order are involved, e.g. vortex lines of length L, then one faces the 
problem that the vortex motion timescale ml? jf? is 5 to 6 orders of magnitude 
too long, when compared to the oscillation period. 

Zero-flow experiments 

Superfluidity is associated with an anomalous mass current response to a gauge 
phase gradient or chemical potential difference. There have been numerous 
attempts in the past to detect superflow in the solid phase with negative results 
(see Ref. [17]). The same conclusion was reached in more recent studies [18], 
both in Vycor and in bulk samples. As they stand now, flow experiments do 
not support, and indeed are incompatible with the supersolid phase of helium. 
We can only note here that in disordered solids large pressure gradients are 
common, and it is not obvious what chemical potential difference is applied to 
the superfluid component in the disordered sample. For example, air pressure 
in the cave deep under the mountain is still at a modest 1 bar value. 



^It should be mentioned that experimental data published so far are rather noisy 
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4.1 Superglass 

Aside looking at defects that are typical for the polycrystalline sample, one 
may wonder if helium can exist in an "ultimate" disordered solid phase, 
namely, a glass. All known structural glasses are in a normal, non-superfiuid 
state. In this regard, ^He offers an intriguing possibility, unique for a quantum 
solid, of being in a "supcrfluid glassy" phase, or superglass (SG). By definition, 
in the metastable SG phase the translation invariance is broken, but the pair 
correlation function g{r) = (n(r)n(O)) features no diagonal long range order. 
At the same time, a superglass has ODLRO and non-zero supcrfluid response. 
Strange as it is, the SG is reminiscent of a sponge soaked in a supcrfluid liquid 
made of the same atoms. 

Superglass was observed in numerical simulations of high-pressure (about 
150 bar) samples prepared by fast temperature quench from the normal liq- 
uid state [32]. Clearly, the dynamics of real helium under cooling is vastly 
different from the thermalization dynamics in the Monte Carlo simulation 
done (i) in imaginary time, and (ii) in the absence of energy conservation. 
Thus, numerical simulations rather answer the question of existence of the 
metastable phase and its quasi-equilibrium properties, but not how easily it 
can be prepared experimentally. One may also wonder if helium SG is merely 
an artifact of the simulation algorithm (though we do not see any obvious rea- 
son to suspect that), or a long-lived, metastable, physical phase responsible 
for the unexpected outcome of the acoustic-pulse experiment [61] which aimed 
at nucleating the crystal phase in the middle of the cold, 50 mK, supcrfluid 
liquid but failed. The local pressure in the pulse reached 160 bar, well above 
the theoretically predicted threshold of 60 bar for quantum nucleation [62] . It 
is not known what the normal-superfiuid transition temperature in the glassy 
phase is, at what temperature the metastable liquid freezes into a glass, and 
what role (if any) the SG might play in the torsion oscillator experiments, e.g., 
by being trapped in some amount between micro-crystals. 

5 Grain boundary superfluidity 

Even if the bulk phase of the material is an insulating solid, there is no physical 
reason why grain boundaries (GB) in the same material ought to be insulating 
as well, especially in the vicinity of the (weakly first-order) liquid-solid tran- 
sition. It means that SF in lower dimensionality can be obtained by simply 
placing two solid pieces right next to each other. Conceptually, this possibility 
was demonstrated in model simulations of domain walls in the checkerboard 
solid formed by hard core bosons with the nearest-neighbor repulsion on a sim- 
ple cubic lattice [63] . Subsequent path-integral Monte Carlo simulations of '^He 
polycrystals [46] did reveal that some (not all !) grain boundaries are likely 
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Figure 4. Phase coherence properties of grain boundaries and ridges in the *He sample consisting 
of two microcrystals with about 12 X 12 X 7 atoms each. Periodic boundary conditions are assumed. 
The XY boundaries between different crystals were created by rotating both of them at random. 

Other grain boundaries are induced by periodic boundary conditions in the simulation box of 
rectangular shape, i.e. they are created within the same crystal. We forced YZ boundaries to be 
insulating by fixing atom positions in the YZ-layer (1.5a thick) at the equilibrium hep points. 
Suporfiudity was then possible only along the x and z directions. Points in the figure are show 
positions of particles which participated in macroscopic exchange cycles with non-zero winding 
numbers, i.e. directly contributed to the superfiuid response. The two panels (from left to right) 
projections on the XY (only the upper half of the sample is projected on the XY plane) and XZ 
planes. The Monte Carlo simulation was performed in the grand canonical ensemble at the 

coexistence curve. 



to support superfluidity. Superfluidity in the layer at the disordered Vycor 
substrate was also reported in Ref. [64]. 

Simulations also find that nearly all ridges (lines of contact between different 
GB) show robust phase coherence properties. It is thus possible that super- 
fluidity across grain boundaries of small (216-atom) crystals is nothing but 
the proximity effect. This issue was addressed in simulations of much larger 
GB [46] in a system of about 2000 atoms consisting of two randomly oriented 
crystallites at the melting density n = 0.0287 . The main conclusions 
did not change; while grain boundaries at special angles and high symme- 
try directions are insulating, generic grain boundaries appear to be superfiuid 
with typical transition temperatures (orientation dependent) of about half a 
Kelvin and the maximum possible Tc at about 1.5 K. The width of the super- 
fluid GB region is about ~ 3a where a is the nearest-neighbor distance in the 
hep crystal, see Fig. 4. 

In a remarkable recent experiment [19], superflow along the grain boundary 
was detected in an experiment whose setup is similar to that shown in Fig. 1. 
The difference between the solid levels was found to decay according to the 
linear law h = const characteristic of the superfiuid flow at the critical velocity 
(estimated to be of the order of meters per second). The flow was detected only 
in the presence of grain boundaries; otherwise the sample did not relax. On one 
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occasion, the dynamics abruptly came to a halt, as the boundary had suddenly 
disappeared. A grain boundary is a topological crystalline defect, and its quick 
(in less than 1 min time) disappearance is by itself an amazing observation, 
since it has no obvious explanation other than suggesting that boundaries are 
in a rough superfluid phase, supporting coherent recrystallization waves. 



6 Perspectives for future theoretical and experimental work 

It is currently our strong belief, that all known experimental facts, together 
with results of first-principles numerical simulations, rule out any explanation 
of the observed supersolid phenomenon in ^He within a homogeneous crystal 
framework. The only reasonable alternative, consistent with many observa- 
tions, is that of supersolidity induced by crystalline disorder. At the moment, 
very little is known for sure about structure properties of solid samples in 
the torsion oscillator experiments. Hopefully, in future investigations it will be 
possible to control sample quality in situ, e.g. by forming solids under con- 
stant pressure, rather than under constant volume conditions. Helium crystals 
of high quality are usually made out of the superfluid liquid at constant pres- 
sure. Solids grown from normal liquid under constant volume conditions are 
not transparent (similar in appearance to milk) [65]. It is also desirable to 
have direct information on the solid order in the same supersolid sample, e.g., 
by optical means. Even more urgently, several experimental observations of 
different groups have to be reconciled, especially regarding pressure depen- 
dence and annealing effects. It seems also that additional measurements of the 
specific heat and other thermodynamic properties can help in quantifying the 
amount of "disorder" in experimental samples [52]. 

Analytic theories of the SFS state in continuous space are phenomenological 
in nature when it comes to predicting, for a particular Hamiltonian, whether 
this state exists and what are crucial parameters for the SFS state in terms 
of bare couplings. Supersolids arise from competition between the kinetic and 
potential energy terms invalidating perturbative and mean-field treatments. 
At present, we do not know a single example of a system characterized by 
a realistic interatomic potential, which has a stable SFS phase in continuous 
space. This work is worth pursuing, since there are exciting new possibilities 
in engineering interparticle interactions in cold atomic systems. It would be 
also interesting to see, once the ideal SFS phase is found, which of the recent 
theoretical predictions capture the microscopic picture of the supersolid phase 
more adequately. 

So far, the only observation of the superfluid glass phase was numerical (at 
T = 0.2 K and n = 0.0359A^'^). Most properties of the helium glass phase, 
and even its very existence, are a mystery. Theoretically, more work has to be 



February 6, 2008 22:41 Advances in Physics perspective 



21 



P 



^ NORMAL GLASS 








s 

IPPRPI 1 lin f^l AQc? ^ * 


NORMAL LIQUID 









Y 



2 T(K) 

Figure 5. A sketch of the metastable hehum glass phase diagram. Apart from the melting pressure 
Pm we do not show the equilibrium liquid-solid line in the figure; the metastable glass state is 
supposed to be inside the stable solid phase. Since the liquid-solid transition is first-order the 
normal-superfluid transition in the metastable liquid is a continuation of the A-line. The dashed 
line is indicating a sharp crossover, not the genuine phase transition; above this line the structural 
relaxation time is beyond the experimental reach. 

done to quantify the SG state, and to see whether the heuristic phase diagram 
shown in Fig. 5 takes place. Experimentally, one may think of preparing glassy 
solid samples by fast density increase at low temperature rather then temper- 
ature quench at high density. Experimentally, it seems that in helium pressure 
increase can be achieved on a much shorter time scale than cooling down to 
T = 0.2 K. [One objection to this proposal is that crystal growth from the 
cell wall at low temperature is extremely fast [66].] Another way of looking 
for non-ergodic properties of helium samples, i.e. system dependence on the 
intial conditions and preparation protocol, is to compare samples cooled while 
rotating with samples rotated after cooling [67]. 

At the moment, the grain boundary flow experiments give us the most un- 
ambiguous clue to the solution of the supersolid puzzle. On the other hand, 
a whole plethora of new questions arise in this connection. The problem is 
very broad, not only because crystalline interfaces can be prepared at dif- 
ferent angles with respect to the crystal axes, but also because one can get 
different results, depending on how well the initial state was "annealed" (with 
respect to minimizing the classical deformation energy), and whether defects 
are "doped" with out-of-equilibrium vacancies and/or intcrstitials (one may 
wonder if vacancies cluster on defects the same way they do it in the bulk and 
whether they can make, say a dislocation core, superfluid ?). 
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In the experimental setup of Ref. [19], the superflow in the boundary layer 
was observed only at the coexistence curve. It is important to extend this 
type of studies to higher pressures, though the experimental design has to be 
modified. One possibility [68], is to look at the same topology as in Fig. 1 
with the bulk solid being sandwiched between two Vycor plates. Since helium 
in Vycor remains superfiuid at pressures about 10 bar higher than the bulk 
melting value Pm one may think of measuring the flow through the solid in 
response to the pressure difference between the two arms of the U-tubc. 

Another fascinating subject is recrystallization waves in the boundary layer. 
Though possible in principle [63] and being a logical explanation for the fast 
disappearance of grain boundaries in the experiment [19], recrystallization 
waves still remain to be seen and detected. Formally, there is no one-to-one 
connection between recrystallization waves and superfluidity in the boundary, 
e.g. smooth domain walls in the 3D checkerboard solid remain superfiuid at 
zero-temperature, but the theory of the interplay between the two phenomena 
is basically non-existent. 

Previous NMR studies of ^He solid solutions [54-56, 59] provided detailed 
information on the isotopic impurity hopping dynamics, as well as on interpar- 
ticle interactions. These studies were performed at relatively high temperatures 
T > 0.7 K and in high quality crystals. It would be interesting to see if NMR 
techniques can be used to understand what happens to ^He atoms in disor- 
dered samples grown from the normal high-density liquid and cooled down 
to the 100 mK range. If ^He atoms form a dense solution at the defects, one 
should be able to see it in the increased NMR relaxation rates since coupling 
between the ^He magnetic moments is the only available mechanism. Another 
exciting possibility is to use NMR techniques to see if ^He solids can be in the 
metastable glass phase. 

The anomalously low critical velocity in the torsional oscillator studies is yet 
to be understood, even at the expense of admitting that the observed phase 
is, in fact, not superfiuid [40]. 



7 Concluding remarks 

Supersolids may happen in two ways: either as ideal incommensurate crys- 
tals with zero-point vacancies or commensurate crystals full of topological 
defects. Though the original observation of supersolidity in ^Hc can not be 
interpreted in terms of zero-point vacancy mechanism in a perfect crystal, 
things turn out to be far more interesting and exciting than expected. Instead 
of one well-characterized crystal, we now face a whole variety of disordered 
solid samples. This immediately brings into focus grain boundaries, disloca- 
tions, and glassy phases which are, on their own right, new superfiuid physical 
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systems. For example, dislocations and grain boundaries should be considered 
as special low-dimensional objects inseparable from their three-dimensional 
environment. Studies of their fascinating properties make the experimental 
and theoretical landscape multi-dimensional, since in addition to pressure and 
temperature, the set of important parameters also includes various orientation 
angles. There is little doubt that this research field is going to last. We also be- 
lieve that superfiuid amorphous helium is a unique and remarkable state which 
will add a new angle to studies of glassy systems. Apart from understanding 
the role of crystalline defects in the supersolid phenomenon, it is extremely 
important to formulate what type of realistic interaction potential between 
particles can lead to the ideal supersolid phase — this long-term goal which 
started more than 50 years ago is yet to be achieved. Existing experimental 
and theoretical tools are more than capable of making significant progress in 
achieving these goals, both at the phenomenological and microscopic levels. 
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